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ABSTRACT 

One of the major sources of the X-ray emitting hot gas around galaxies is the feedback from 
supernovae (SNe), but most of this metal-enriched feedback material is often not directly detected 
in X-ray observations. This missing galactic feedback problem is extremely prominent in early-type 
galaxy bulges where there is little cool gas to make the Type la SNe ejecta radiate at lower temperature 
beyond the X-ray domain. We herein present a deep Suzaku observation of an SO galaxy NGC 5866, 
which is relatively rich in molecular gas as an SO galaxy and shows significant evidence of cool-hot 
gas interaction. By jointly analyzing the Suzaku and an archival Chandra data, we measure the Fe/0 
abundance ratio to be T.SSlg 5 ® relative to solar values. This abundance ratio is much higher than 
those of spiral galaxies, and even among the highest ones of SO and elliptical galaxies. NGC 5866 
also simultaneously has the highest Fe/0 abundance ratio and molecular gas mass among a small 
sample of gas-poor early-type galaxies. An estimation of the Fe budget indicates that NGC 5866 
could preserve a larger than usual fraction, but far from the total amount of Fe injected by Type la 
SNe. We also find that the hot gas temperature increases from inner to outer halos, with the inner 
region has a temperature of ~ 0.25 keV, clearly lower than that expected from Type la SNe heating. 

This low temperature could be most naturally explained by additional cooling processes related to 
the cool-hot gas interaction as being indicated by the existence of many extraplanar dusty filaments. 

Our results indicate that the large cool gas content and the presence of cool-hot gas interaction in 
the inner region of NGC 5866 have significantly reduced the specific energy of the SN ejecta and so 
the velocity of the galactic outflow. The galaxy could thus preserve a considerable fraction of the 
metal-enriched feedback material from being blown out. 


1. INTRODUCTION 

The hot X-ray-emitting coronae around disk galaxies serve as a reservoir where galaxies acquire baryons to build up 
galactic disks and deposit kinetic ener gy and chemically-enriched materials. This X-ray emitting gas has been detected 
in various types of disk galaxies le.g.. lLi et al.ll200M I20IH [Li fc Wan Bll2013^ iMineo et al.l|201 ^ . and can be globalR 
reproduced by recent cosmological hydrodynamical simulations ( Crain et al T l2mnHLi et al.ll2m4HBogdau et al.llWl 
supp o rting the expected origin of galactic corona from feedback close to the galaxy (jStrickland et al.l 12004 iLi et al.l 
20M iLi fc Wandboi.lbH. or maybe related to accretion of ext ernal gas at large galactocentric radii (r > 50 kpc; 


Anderson fc Bregman Il201lt iDai et al1l2012l : iBogdan et al.l 1201^ . 


The coronal properties are observed to correlate with various galaxy properties l|Li fc Wang|l2013bl) . In particular, 
there is a tight correlation between the coronal soft X-ray luminosity and the total [core collapsed (CC)-l-Type la] 
supernova (SN) energy injection rate. The X-ray radiation efficiency (the fraction of SN energy observed in the coronal 
X-ray emissio n) is quite low (^0 .4%), indicating that most of the feedback material is missing from a clear detection 
in soft X-ray (ILi fc Wang|l20I3bll . In galaxies rich in cool gas, a large fraction of the feedback ene rgy may be radiate d 
through the recombination line of warm ionized gas at a lower temperature {kT ~ 10^“® K; e.g.. lAgertz et al]l20I,‘lfl . 
However, in early-type galaxy bulges which are typically poor in cool gas, the missing feedback mate rial is most likely 
carried out to the intergala.ctic sp ace by a galactic bulge wind mainly produced by Type la SNe ('e.g.. lTang et al.ll2009l : 
iTang fc Wang|[2M^ lWang|[2MCT . 

The metal-enriched galactic feedback material is typically hot {kT > 10^ K) and has t oo low density (n < 
a few X 10“^ cm“^) to be directly detected in current X-ray observations (|Tang et al.l 120091 : iTang fc Wand [20T(»1 . 
It is thus controversial what is the nature of the hot gas detected in soft X-ray around near by galaxies: the real di ffuse 
hot tenuous corona or just the mixed material at the interface of cool and hot gases le.g.. iStrickland et al.ll20(i^ ? In 
order to understand the nature of this soft X-ray emission and to search for the missing galactic feedback material, 
we need to measure the metallicity of hot gas around isolated, moderate- or low-mass, early-type galaxies, which are 
weakly affected by the intra-cluster medium (ICM), gravitational heating, or current star formation. However, existing 
X-ray observations of such galaxie s often reveal a low Fe (a primary product of Type la SN) abundance, typically 
subsolar or slightly supersolar ([e.g.. lHumDhrev fc Buotel[2006l : lJi et al.ll2009l : iKonami et al.l[^0I4ll . indicating that most 
of the SN ejecta is not detected. 

We herein present a deep Suzaku observation of NGC 5866. This galaxy locates in a low density environment, with 
two companions, NGC 5907 and NGC 5879, located several hundred kpc away in projection; their interactions with 
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NGC 5866 should have negligible effects on the coronal properties. The AGN and star formation activities are both 
weak, unlikely contribnting significantly to the X-ray emission. The edge-on orientation helps the gas-rich galactic disk 
to block most of the on-disk X-ray-bright point-like sources, reducing the complexity of analyzing th e low-resolutio n 
Suzaku data. Based on our preliminary study of NGG 5866 using Chandra and multi-wavelength data (jLi et al.ll2009fl . 
we found it has high Fe/0 abundance ratio as a disk galaxy (iLi &: Wang|l2013blf . thus well-suited for studying Type la 
SN metal enrichment. However, the metal abundances are poorly constrained with Chandra data only, so we proposed 
deep Suzaku observation of this galaxy to collect more photons 100 ks; Gycle 6; PI: Li J. T.). 

In this paper, we jointly analyze the Suzaku and Chandra data, and further compare this galaxy to other galaxies 
and discuss how the new observation could help us to understand the Type la SN feedback processes. The paper is 
organized as follows: In Sj2l we describe the calibration of the Suzaku data and the joint analysis of the Suzaku and 
Chandra spectra. In ^ we compare our X-ray measurements of NGC 5866 to similar measurements of other galaxies 
and further discuss the results and their implications. Out major results and conclusions are summarized in 21 Errors 
are quoted at 90% confidence interval throughout the paper, unless specifically noted. 

2. DATA REDUCTION AND ANALYSIS 
2 .1. Suzaku data reduction 

Suzaku observation of NGC 5866 was taken on 2011-05-20 with a total exposure time of 101.8 ks. We reprocess the 
XIS data using the ftool aepipeline and CALDB files of version 2014 October 1. We identify the high background 
flare time intervals using a 3 cr clip with a lightcurve binned at 200 s. We then remove these background flares and 
obtain an effective exposure time of 100.8, 80.9, 100.2 ks for XISO, XISl, XIS3, respectively. We generate the XIS 
non-X-ray background spectra and images using the ftool xisnxbgen. In imaging analysis, the vignetting is corrected 
by simulating an event file using xissim. The images are rescaled with bin=8 (pixel size=8.3"), trimmed to remove 
low exposure regions, flat fielded, and exposure corrected. The resultant tri-color images in 0.3-1 keV, 1-2 keV, and 
2-7 keV are shown in Fig.[T^, while the 0.3-1 keV contours in the central region are presented in Fig.[Tl3. 
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Fig. 1.— (a) Suzaku tri-color (red: 0.3-1 keV; green: 1-2 keV; blue: 2-7 keV) counts image of NGC 5866 and the surrounding regions. The 
annulus marks the region used to extract the sky background spectrum, after removing the sources in the two smaller circles. The small box 
marks the central 4' x 4' region centered at the nucleus of NGC 5866, and is also the FOV of panel (b). (b) Suzaku background-subtracted 
and exposure-corrected 0.3-1 keV contours overlaid on the SDSS dr7 tri-color images (red: i-band; green: r-band; blue: g-band). The inner 
circle has a radius of 1' and is used to extracted the inner halo spectrum, while the outer circle has a radius of 1.6' and the annulus between 
these two circles is used to extract the outer halo spectrum . The small circles mark the location and size (96% PSF energy enclosed radius) 
of the Chandra detected point-like sources lILi et al.ll2009l) . 


Suzaku XRT has a typical size of the point spread function (PSF) of ~ 2'. As shown in Fig. [T}a,, except for the 
emission from NGC 5866, the other bright X-ray sources (the two sources at the edge of the CCD and the two fainter 
sources marked with the circles) are all background sources also identihed by the Chandra ( W2.2\ iLi et al.ll2009f) . 

We define the regions used for extracting sky background and source spectra as shown in Fig. [IJi,b. The smaller 
circle in Fig.[T|D has a diameter of 2', matching the XRT PSF size; we herein call it the “inner halo” region. The larger 
circle in Fig. [TJd has a diameter of 3.2'; we define the annulus between these two circles as the “outer halo” region, 
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TABLE 1 

Hot gas properties from spectral analysis 


Parameter 
X^/d-o.f. 
kT (keV) 

O abundance (solar) 
Fe abundance (solar) 
EM (cm-^^kpc^) 


Inner halo 


121 . 00/135 

0 . 2481 °;“? 

0.2ll°;l° 

1 r7+l-32 

-0.58 

0 . 053+°“5 


Outer halo 


6 . 59/15 


0.015 ± 0.004 


Spectral analysis results of the “inner halo” region is obtained by jointly fitting the Suzaku and Chandra spectra (in total 10 free parameters). 
Spectral analysis results of the “outer halo” region is obtained by fitting the Suzaku XISl spectrum with a model in which the absorption 
column density and photon index of the power law, as well as the O and Fe abundances of the hot gas, are all fixed at the values of the 
inner halo (in total 5 free parameters). See text for the detail description of the spectral models and Fig. [2^,c for the fitted spectra. All 
the errors quoted in this table are statistical only, at 90% confidence level. EM is the emission measure of the VAPEC component. 

outside which there is no significant diffuse X-ray emission. The large annulus in Fig. [TJi (with diameters of 6' — 10'), 
after removing the two bright background sources, are used to extract the sky background spectrum. 

We extract the XIS spectra and the corresponding response files of the source and sky background regions using the 
ftool xselect. The non-X-ray instrumental background spectra are also extracted from the same regions. There is a 
weak line feature at ~ 0.62 keV of the XISl spectra, which is not detected in XIS0/XIS3 or the Chandra spectra. This 
feature is not common in thermal plasma spectrum. We do not find any significant time variation of it. The feature 
also presents over the entire XISl field of view. It is thus a background feature most likely caused by quiescent solar 
wind charge exchange (e.g., lYoshitake et al.ll2013D . We will use an additional Gauss model to account for it in the 
following spectral analysis. 


2 . 2 . Archival Chandra data 

Suzaku has too low angular resolution to individually detect X-ray bright point-like sources, such as AGNs or X-ray 
binaries, within the optical extent of NGC 5866. Unlike X-ray faint stellar sources, the contribution from these X-ray 
bright sources is quite uncertain, and could be the major uncertainty in characterizing the truly hot gas emission. We 
therefore need a joint analysis of the Chandra and Suzaku data to better decompose the X-ray emission from hot gas 
and stellar sources and/or AGN. 

The archival Chandra ACIS-S observation of NGC 5866 is presented in ILi et HI (12009112011[l : iLi &: Wand l)2013aD . 
It has an effective exposure time of ^ 24.3 ks. We herein directly use the calibrated data, the diffuse X-ray images, 
and the detected point-like sources (Fig. [Ho) from ILi k. Wand (|2013al) . The Chandra spectra are extracted from the 
same regions as the Suzaku spectra. We do not remove the point sources from the diffuse X-ray emission, in order to 
be consistent with the Suzaku spectra. Instead, we extract an accumulated Chandra spectrum of all the point sources 
detected within the inner halo region (point sources in the outer halo have too few counts to construct a spectrum). 
This spectrum will be jointly analyzed with the Chandra/Suzaku diffuse X-ray spectra extracted from the inner halo. 

2 . 3 . Spectral analysis 

In general, the X-ray spectra of the inner and outer halos are comprised of several different components. We 
characterize the hot gas emission with an optically thin ionization equilibrium thermal plasma model VAPEC. The 
oxygen and iron abundances are set free, while the abundances of other elements are set to the solar values. In 
order to better characterize the hot gas emission, we also need to quantify the various stellar contributions. With the 
archival Chandra data, we reach a point source detection limit of ~ 4 x 10^^ ergs s“° in 0.5-8 keV. The accumulated 
spectrum of all the detected point sources within the inner halo region can be well modeled with a power law, which 
has a photon index of F = 1.4 and an absorbing column density slightly larger than the Galactic foreground value 
of Vh ^ 1.46 X 10^° cm“^. Below the detection limit, the unresolved stellar source emission is dominated by the 
emission from cataclysmic variables (CVs) and coronal active binaries (ABs), as well as the emission from unresolved 
lo w mass X - ray b inaries (LMXBs) and the residual photons spilling out of the discrete source removal circles. Similar 
as ILi et al.l (1200^ . we model this component with a 0.5 keV MEKAL model (optically thin thermal plasma) and a 
F = 1 .9 power law model. B oth components are well calibrated and their normalizations are scaled with the stellar 
mass (jR.evnivtsev et al]l2008f l. There is still a small residual at ~ 0.62 keV after subtracting the sky background. 
This feature is only presented in the XISl spectrum and is most likely a residual feature produced by the quiescent 
solar wind charge exchange. We therefore add into the XISl spectrum a narrow gauss model with fixed line width of 
10“^ keV to account for it. 

The Suzaku (XISO, XISl, XIS3) and Chandra (total emission and detected point source emission) spectra are jointly 
fitted, with most of the model parameters fixed. Fitting the absorption column density (Ah) of the thermal component 
results in an even lower Ah than the foreground value and comparable values of other parameters and x^/d.o./.. We 
therefore fix it at the foreground value of Vh 1-46 x cm“^. The only free parameters are the temperature, O 
and Fe abundances, normalization (or flux) of the hot gas, the absorption column density, power law photon index and 
normalization of the Chandra-detected point sources, as well as the centroid energy and normalization of the gaussian 
line added to the XISl spectrum. The joint analysis gives us much better counting statistic than presented in previous 





























4 




Fig. 2. — (a) Joint analysis of the Suzaku XISO (black), XISl (red), XIS3 (green), Chandra total spectrum (blue), and Chandra spectrum 
of point-like sources (cyan) from the “inner halo” region. Different model components of the XISO spectrum, which has the best counting 
statistic, are also plotted as red dash-dotted curves, (b) Confidence range of the O and Fe abundances of the model in panel (a) [contours 
at 68.3% (1 fj), 90%, 95%, and 99% confidence levels]. The plus sign marks the best-fit values. The dashed lines mark different values of 
Fe/O abundance ratios, (c) XISO spectrum of the “outer halo” region, with different model components marked with dash-dotted curves. 
The abundances of O and Fe, and the power law photon index are all fixed at the values of the “inner halo”. 


works (|Li et al.ll2009L [2M1I : iLi fc Wandl2013aD . The best-fitted spectral model parameters are summarized in Table [T] 
and the fitted spectra are presented in Fig. [2^. One of our primary goal in this paper is to better determine the Fe/O 
abundance ratio. We plot the confidence range of Fe and O abundances in Fig. [2b to further examine how reliable 
is our measurement. It is highly reliable (at > 90% confidence level) that the hot coronal gas in the inner halo of 
NGC 5866 has a super-solar Fe/O abundance ratio (Fe/O = 
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Fig. 3. — This figure shows how the assumptions on the power law component and the different values of Fe/O ratio may affect the 
measurement of the hot gas temperature {kT) in the “outer halo”. The black, green, and blue symbols are for models 1-3 as described in 
the text. They are at the same Fe/O value marked by the black symbols, but the green and blue symbols have been shifted a little along 
the Fe/O axis to clarify. The red dot and error bars show the measurement and 90% confidence error of the “inner halo” region. 


For the outer halo, only the XISl spectrum has enough counts for spectral analysis. We fit this spectrum with the 
same model as for the inner halo. However, the O and Fe abundances of the hot gas are poorly constrained. We 
then fix the abundances at the values of the inner halo, in order to estimate the hot gas temperature. The resultant 
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X^/d.o.f. is too low (~ 0.44; Table [T]), indicating that there may be too many free parameters (five, temperature and 
normalization of the thermal plasma, normalization of the power law component, as well as the centroid energy and 
normalization of the gaussian line) and large systematic uncertainties. 

The largest systematic uncertainty which could highly affect the characterization of the hot gas emission may be the 
presence of the poorly constrained power law component. As the most important parameter which can be determined 
in the outer halo is the hot gas temperature, we further examine how the assumption of the power law model may affect 
the measurement of the temperature. We fit the spectrum in three different ways: (1) the absorption column density 
and photon index of the power law are fixed at the values of the inner halo; (2) set a lower limit on the absorption 
column density (at the Galactic foreground value), but free the photon index; (3) remove the power law component 
and fit the spectrum only in 0.3-1.5 keV, where the X-ray emission is dominated by the thermal component. In each 
model, we also test how the results are affected by different assumed values of the Fe/0 abundance ratio (change 
the Fe abundance while fix the O abundance at 0.2 solar). The results are shown in Fig. [U in comparison with the 
temperature of the inner halo. It is obvious that in any models, the hot gas temperature of the outer halo is higher 
than that of the inner halo. 

2.4. Systematic uncertainties of the spectral analysis for the inner halo 

The errors quoted in Table [1] only include the statistical uncertainties of spectral fitting. We herein discuss some 
systematic uncertainties of the temperature and abundance measurements in the inner halo. 

The contribution of the CV-I-AB component can exhibit a factor of ^ 1.7 vari ation in different galaxi es (e.g., 
iBoroson et al.l [2?iTl^ . We have adopted a relatively high Tx,cv-i-AB/Af* ratio from iR.evnivtsev et al.l (j2008ll . which 
may affect the measured temperature and metallicity. If we adopt a CV-I-AB model which has only ~ 50% of the 
flux of the current model, the measured temperature and metallicity are only slightly different: kT = 0.255^q q 27 ! O 
abundance= 0.20to,o9, Fe abundance= 1.44to;||. Therefore, the small uncertainty in the CV-I-AB component could 
not significantly affect our results. 

We have assumed a single -temperature model of the therm al plasma, which is in general true for early-type disk 
galaxies like NGC 5866 le.g.. lLi et al.ll20lII : ITT fc WanBll2013a[l . However, it is still possible that the hot gas has a more 
complicated thermal structure. The low best-fit temperature of the hot gas indicates that it is most likely produced 
in the interface between the cool gas and a volume-filled hot tenuous galactic outflow (see the discussions in 113.211 . We 
further discuss how the presence of a high-temperature component representing the galactic outflow may affect our 
measurements of the metallicity of the low-temperature hot gas. The properties of this galactic outflow component is 
poorly known observationally. We therefore adopt the gas properties i nferred from numerical simulations of a galactic 
bulge without much cool gas (|Tang et al.l 120091 : iTang fc Wand [20T(ill . According to these simulations, the hot gas 
temperature has a broad distribution, but the emission measure or radiative cooling coefficients typically peak at 
^ 1 keV. We add a VAPEC model with a fixed temperature of 1 keV to the model of the inner halo. The quality 
of the spectra do not allow us to measure the metallicity of the two components separately. We therefore simply link 
the O and Fe abundances to those of the low-temperature component. The best-fit model indicates an even lower 
temperature of the low-temperature component (^ 0.18 keV), and an even higher Fe metallicity of > 10 solar, further 
strengthening our arguments for the low-temperature and high-Fe-abundance of the hot gas in the inner region of 
NGC 5866. 

The applied index of the power law component (~ 1.4), as determined with the Chandra resolved point sources, 
is probably too low for the LMXBs in the galact ic bulge. In nearb y early-type galaxies, this index typically span in 
a range of 1.4-1.8, with a best fit value of 1.56 l)Irwin et al.l 120031) . We test how the value of this power law index 
may affect our results by fixing it at 1.8. The resultant best-fit temperature and metallicity of the inner region are: 
kT = 0.232^0.033, O abundance= 0.20lo |!)g, and Fe abundance= consistent with the results obtained by 

directly fitting this power law index. 

Based on the discussions on the systematic uncertainties of the spectral analysis for the inner halo, we conclude that 
NGC 5866 has a significantly super-solar Fe/0 abundance ratio and a low hot gas temperature of < 0.25 keV in the 
inner region. 

3. DISCUSSION 

3.1. Metallicity of the hot halo gas 

iLi &: Wangl (|2013bll found an anti-correlation between the hot gas Fe/0 abundance ratio and the morphological type 
code (TC) of the galaxy, indicating that Type la SN has a larger contribution in metal enrichment in earlier-type 
galaxies (smaller TC) . We herein compar e the hot gas metallicity of NGC 5866 to an updated version of the TC-Fe/0 
diagram presented in ILi fc Wand ()2013bf l. As shown in Fig. |4j NGC 5866 is clearly richer in Fe than other spiral and 
SO (TC ^ —3 — 0/ galaxies, and even many of the earlier-type elliptical galaxies. Since most of the metal-enriched 
mate rial in gas-poor SO or elliptical galaxies is thought t o be carried out by a galactic scale bulge wind or outflow 
le.g.. lTang et al.ll2009l:[Tang fc Wangjl201(iILi et aT1l2011f) . this high Fe abundance means there is a larger than usual 
fraction of Type la SN ejecta detected in soft X-ray in the inner halo of NGC 5866. 

We estimate the total mass of Fe detected in NGC 5866. Assuming spherical symmetry of the X-ray emitting regions, 
we could estimate the density of hot gas from the spectral fitting parameters. The densities of the inner and outer halos 
are ~ 1.2 x 10“^/“^/^ cm“^ and ^ 3.6 x 10“^/“^/^ cm“^, respectively, where / is the poorly known volume filling 
factor. The total mass of hot gas contained in the inner and outer halos are both ^ 10®/^/^ Mq. As the metallicity 
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Fig. 4.— Fe/O abundance ra tio v.s. the morphological t ype code (TC) of NG C 5866 (black hou rglass), compared to the X -ray measure- 
ments from [Li &: Wan^ Il2013al[b|) (red and black symbols). IMartin et al.l (120021) (green diamond). IHumphrev fc Buotd (120061) and IJi et al.l 
([ 2003 ) (blue symbols). The error bars of the Fe/O ratio are plotted at 1 cr confidence level to clarify. 

in the outer halo cannot be well constrained, we only estimate the total Fe mass contained in the inner halo, which is 
- 2 X 10^ M0. 

We further estimate the expected metal (Fe) budget in NGC 5866. In isolated gas-poor elliptical or SO galaxies, 
the major source of gas is from stellar mass loss while the major sourc e of metal is from Type la SNe. Adopting a 
Type la SN rate of 4.4 x 10“^ SN yr“^/10^° Mq (jMannucci et al.l|2005ll and the stellar mass of 4.5 x 10^° Mq in the 
inner halo (estimated from the K-band luminosity), we could obtain a total Type la SN rate of 2 x 10“^ SN yr“^. 
If no SN ejecta can be confined around the galactic bulge, i.e., they are carried out by a galactic outflow, we would 
expect the detected metal is only produced within the dynamical cross timescale of the outflow. Assuming a typica l 
outflow velocity of ^ 300 km s“^ from numerical simulations of a similar gas-poor galactic bulge (|Tang et al.l 120091) . 
the radius of 1' of the inner halo gives a dynamical cross timescale of tdyn ^ 1-5 x 10^ yr. Within this timescale, 
there are 2 x 10^ Mq Fe injected into the ISM by Type la SNe (assuming each Type la SN eject ~ 0.7 Mq Fe; 
iNomoto et al.l [198411 . about one order of magnitude smaller than what we have detected in the inner halo, assuming 
the volume Filing factor / 1. On the other hand, if most of the SN ejecta can be confined around NGC 5866, the 

total amount of metal injected by SNe should be calculated within a timescale since the last major starburst stage. 
This timescale can be estima ted with t he glo bular cluster (GC) specific frequency. We herein adopt the age after last 
starburst of NGC 5866 from iLi et al.l (|2009tl : tage 8 Gyr. Within this timescale, the total Fe injected by Type la 
SNe is 4^ 10^ Mq, about 50 times higher than what we have detected. 

The most nncertain parameter in the above estimation is the volnme filling factor /. In starburst galaxies, the volume 
filling factor is often extremely low (e.g., < 2%) due t o the strong interaction bet ween the fast galactic superwind and 
the large amount of highly-structured cool gas le.g.. iStrickland fc Stevendl200fl[ l. In star formation inactive galaxies 
such as NGC 5866, it is unlikely that / is such extreme, because the amount of cool gas is much less and there is much 
weaker cool-hot gas interaction. We do not have a direct observational constraint on /. However, even if we adopt 
the extreme volume filling factor of starburst galaxies, the estimated total Fe mass is still consistent with or slightly 
larger than the Type la SNe products within the dynamical cross timescale of the outflow. Therefore, we conclude 
that Type la SNe should be sufficient to explain the observed metallicity of the hot gas in NGC 5866. The relatively 
larger Fe abundance of NGC 5866 compared to other elliptical or SO galaxies indicate that NGC 5866 is able to keep a 
larger fraction of the SNe ejecta surrounding the galaxy, although most of the metals are still blown out, as indicated 
by the large difference between the observed and expected metals during the lifetime of the galaxy. 

3.2. The effect of cool gas 

The perfect edge-on orientation of NGC 5866 allows us to search for dusty filaments through absorption features 
on th e optical images, whic h is the best evidence for the presence of extraplanar cool gas in such gas-poor galaxies 
('e.g.. lHowk fc Savagelll999ll . We present high resolution unsharp-masked HST images of NGC 5866 in Fig. [51 which 
show prominent dusty filaments extended perpendicular to the galactic dis k. The largest filaments extend as far as 
4^ 20" ( 4 ^ 1.5 kpc) above the disk (marked by the white arrow in Fig. |5|). iLi et al.l (|2009ll found that most of these 
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dusty filaments, especially those close to the galactic disk, could be produced by single SN. Therefore, Type la SN is 
energetically sujficient to power the cool-hot gas interaction in this SF inactive SO galaxy. 
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Fig. 5.— The unsharp-masked HST tri-color images (red: F625W; green: F555W; blue: F435W) of the central 1.2^ x 0.8^ of NGC 5866, 
showing the dusty (dark) extraplanar filaments. The arrow marks one of the largest dusty filaments, which is ~ 20^^ in length perpendicular 
to the galactic disk. 


We further investigate the role of cool gas on the metal enrichment by comparing the Fe/0 abundance ratio and the 
total molecular gas mass (Mhj) of a small sample of gas-poor early-type galaxies (Tabled]). Most of these galaxies 
are too gas-poor so the listed Mhj are upper limits, except for NGC 5866 and NGC 1399. It is therefore difficult to 
argue for a clear positive correlation in Fig. |6l Nevertheless, it is clear that NGC 5866 is the most gas-rich galaxy in 
this sample and shows the highest gas-phase Fe/0 abundance ratio, possibly indicating that its high cold gas content 
may play a role in retaining the metal-enriched SN ejecta. 

The temperature of the hot gas in the inner halo seems too low to be primarily heated by SNe. In gas-poor galactic 
bulges like NGC 5866, hot gas is the d ominant phase an d the radiative cooling of the shocked ISM is overall not 
important in the SN energy budget fe.g.. lLi M. et al.l[20I5[) . The specific enthalpy of old stellar feedback (mass loss of 
evolved low-mass stars and Type la SNe heating) is ~ 4 keV per particle. Such a specific enthalpy gives an average 
gas temperature of ^ 1 keV, considering the energy equipartition bet ween the hot gas thermal energy and the energy 
store d in magnetic field and cosmic rays as well as turbulent motion (|Li et al.l 120091 iTang et all 120091 iTang fc Wangl 
I2010f) . On the other hand, adopting a central velocity dispersion of (t„ ^ 162 km s“^ obtained from the HyperLeda 
database, the expected temperature of gas thermalized in the gravitational potential of NGC 5866 would be ~ 0.2 keV, 
quite close to the measured hot gas temperature of ~ 0.25 keV in the inner halo (Table]!]). Therefore, the low measured 
hot gas temperature apparently indicates that any additional sources of non-gravitational heating, such as Type la 
SNe, play a relatively minor role (contributes only ^ 20% of the total thermal energy). 

Although located in a very low density environment (the environmental gala xy number density of NGC 5866 is 
p = 0.24; in comparison, p > 0.6 is classified as clustered in iLi fc Wang||^13a(l . it is still possible that there exists 
a low-density intra-group medium around NGC 5866. However, neither of the existing Suzaku or Chandra data as 
presented in this paper is deep enough to probe the hot gas distribution at large radii. The Chandra X-ray intensity 
profile s at small radii, appr oximately within the spectral analysis regions of the present paper, are presented in ILi et ahl 
(|2009fl : lLi fc Wan3 (l2013all . These profiles do not reveal any diffuse soft X-ray emission above the stellar contributions 



























TABLE 2 

Hot gas metallicity and molecular gas mass of early-type galaxies 


Name 

Fe/O 

Fe/O (4') 

log Mh 2 /Mq 

Mh 2 ref. 

NGC 5866 

7 63+’^-^^ 

' •'^'^-5.52 

- 

8.47 ±0.01 

Youngll 

NGC 720 


U.1U_1 30 

< 7.38 

Sage07 

NGC 741 

2.70tfH 

- 

< 8.28 

Evans05 

NGC 1399 

9 qi -1-0.31 

1 QO+0.03 

i.OU_o.o4 

7.47 

PrandonilO 

NGC 1407 

2.70lj;«3 

- 

< 7.53 

Sage07 

NGC 4365 

9 1 7 -I-I. 7 O 
< -8.32 

- 

< 7.62 

Youngll 

NGC 4406 


1 qoLU.Ub 

i.OO_Q 05 

< 7.40 

Youngll 

NGC 4472 

2.06l°i® 

2.06t°:°® 

< 7.25 

Youngll 

NGC 4552 

1 4C:;+0.45 

-L-^0_q 43 

3.00t°;« 

< 7.28 

Youngll 

NGC 4636 

-L-D1 _o.i6 

i.59t°;gl 

< 6.87 

Youngll 

NGC 4649 

9 Qr+0.56 

Z.OO_o.60 


< 7.44 

Youngll 

NGC 5846 

5.0011:^® 

- 

< 7.78 

Youngll 


Fe/O abundance ratio of hot gas and molecular gas mass (Mn^) of NGC 5866 and the elliptical galaxies from [Humphrey &: Buot^ (I20Q6I) 
and IJi et al.l l|2QQSIL For galaxies included in bothlHumohrev fc Buotj Il2nnfil^ and IJi et alJ 1|20Q9IL we adopt [Ti et alJ ll2009IV s measurements 
as they use better X-ray data. IJi et al.) (120091) perform X-ray measurements based on different instruments {Chandra ACIS, XMM- 
Newton EPIC and RGS) and with different aperture sizes (diameters of 1' and 4'). We herein adopt the XMM-Newton EPIC P-diameter 
measurements, which should be the most similar as our measurement of NGC 5866, and is less affected by the c l ustere d environment of 
some galaxies. The Fe/O ratio used in the plots are summarized in the “Fe/O” column. We also list IJi et al.l (120091) ’s XMM-Newton 
EPIC 4'-diameter measurements in the “Fe/O (4^)” column for comparison. The differences in these two columns do not affect our results 
discussed in this paper. Errors of the Fe/O ratio are quoted at 90% confidence level in this table, but are plotted at 1 cr confidence level 
in Figs. [4| and [6] to clarify. The molecular gas mass is based on CO observations, with references listed in the last column. Evans05: 

lEvans et al.l l|2005l) : SageOT: ISage et ^ (|2007l ): PrandonilO: IPrandoni et al.l l|201(]|) : Youngll: lYoung et al.l (|2011l) . 



(^solar) 


Fig. 6.— Fe/O abundance ratio v.s. the molecular gas mass {Ma^) of NGC 5866 compared to elliptical galaxies from lHumohrev &: Buot^ 
([20061) and □r et al.l l|2009l) . Except for NGC 5866 and NGC 1399, all the other measurements of are upper limits {Mn^ and the 

references are summarized in Table n. The error bars of Fe/O ratio are plotted at 1 cr confidence level to clarify. 


^ r (or roughly the inner halo region in this paper; Fig. (Hd). Therefore, based on the existing X-ray data, we do 
not find significant evidence for the exis tence of an intra -group medium, which, if presents, could help to raise the hot 
gas temperature in the outer halo ie.g.. lSun et al.l[2?i?)^ . 

It is suggested that the X-ray emissivity of early-type galactic coronae, and also the fraction of SN ejecta detected in 
soft X-ray, a re highly affect ed by the specific energy of the hot gas (energy per unit mass, or temperature if the gas is 
thermalized: iLi et al.ll201ll) . If this specific energy is reduced by some processes, there will be a lower temperature of 
the hot gas and slower galactic outflow. As a result, the hot gas density and X-ray emissivity will be raised and there 
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will be a larger amount of heavy elements produced by Type la SNe (such as Fe) detected in soft X-ray. Cool-hot gas 
interaction, as indicated by the relatively large amount of cool gas (Table [SJ Fig. |6l) and the presence of extraplanar 
dusty hlaments in NGC 5866 (Fig. [Sj, could provide a natural explanation of both the low temperature and high Fe 
abundance of the hot gas. 

The apparently low Type la SNe contribution to the heating of the hot gas is very unlikely, because NGC 5866 shows 
so many features related to Type la SNe feedback, such as the Fe enrichment detected in hot gas and the consistency of 
the gravitational energy of many extraplanar dusty filaments with a single SN energy (|Li et al.ll2009ll . Furthermore, the 
significant raisin g of the hot gas temper ature in the outer halo (Fig. [Sj 112.31) in such a low-mass isolated galaxy is quite 
uncommon fe.g.. iFukazawa et al.ll20d^ . Although we cannot absolutely exclude the contribution from a low-density 
intra-group medium based on the existing Suzaku or Chandra data, it is natural that there exists some additional 
cooling processes, such as the cool-hot gas interaction. The presence of cool gas could affect the surrounding hot 
gas properties via mass loading and enhanced radiative cooling, both processes reduce the specific energy and so the 
temperature of the hot gas. 

The unusually strong cool-hot gas interaction in NGG 5866 as an SO galaxy could also help to explain the high Fe 
abundance detected in soft X-ray. The cool gas is often less metal-enriched, and the mass loading of cool gas into 
hot phase will result in a lower average metallicity of the hot gas. However, most of the SN ejecta contained in the 
high-temperature (typically > 1 keV) hot gas has too low emissivity to be detected in soft X-ray, because the hot 
gas density of gas-poor galactic bulge is too low (typically < 10“^ cm“^, e.g., iTang &: Wang|[2M^ . Therefore, the 
cool-hot gas interaction will have an opposite effect on the observed gas metallicity, by significantly increasing the 
soft X-ray emissivity and so the emission-measure-weighted metallcity as obtained from X-ray observations. Such a 
differ ence between the emission-measure - and mass-weighted metallicity is also evidenced by numerical simulations 
le.g.. iTang et al.ir2009l : iGrain et al.ll2013ll . indicating that we ar e mi ssing from the detection in X-ray observations of 
a large fraction of gas-phase metals, as also being estimated in 113.11 

4. SUMMARY AND CONCLUSION 

We present detail calibration of the Suzaku observation of an isolated SO galaxy NGG 5866. By jointly analyzing this 
Suzaku data together with an archival Chandra observation, we obtain the hot gas temperature in the inner (r < 1') 
and outer (1' < r < 1.6') halos, as well as the O and Fe abundances in the inner halo. These measurements of the 
hot gas properties are much more accurate than previous measurements with the Chandra data only, allowing for a 
comparison with other galaxies. 

We find a significant supersolar Fe/0 abundance ratio of Fe/0 = T.OSj'jg'g® at 90% statistical confidence level in 
the inner halo of NGG 5866, indicating a clear metal enrichment from Type la SNe. The Fe/0 abundance ratio of 
NGG 5866 is also significantly higher than other disk (spiral-|-S0) galaxies, and even many of the elliptical galaxies. 
We also find that NGG 5866 simultaneously has the highest Fe/0 abundance ratio and molecular gas mass {Mh^) 
among a small sample of gas-poor early-type galaxies, although the poor and mostly upper limit measurement of Mhj 
prevents us from claiming a clear positive correlation between Fe/0 and Mhj. The total amount of Fe detected in the 
inner halo of NGG 5866 is 2 x 10®/^/^ Mq. In contrast, the total amount of Fe injected by Type la SNe within the 
dynamical cross timescale of the galactic outflow is 2 x 10^ Mq, while that within the age of the galaxy since the 
last starburst is ^ 10^ Mq. Gonsidering the uncertainty of the volume filling factor / which cannot be too low in this 
non-starburst SO galaxy, we could conclude that NGG 5866 could preserve a larger than usual fraction, but far from 
the total amount of Fe injected by Type la SNe. 

NGG 5866 exhibits a lot of extraplanar dusty filaments extending perpendicular to the galactic disk within the inner 
halo. These filaments indicate the existence of interaction between SNe and cool gas, which probably triggers cool 
gas mass loading into the hot phase and enhances the radiative cooling. We find clearly lower hot gas temperature in 
the inner halo than in the outer halo, which is apparently a result of this cool-hot gas interaction, if this well isolated 
galaxy is not highly affected by a high-temperature intra-group medium. This low temperature in the inner halo is 
also roughly consistent with gravitational heating of the material lost by evolved stars, apparently leaving little room 
for additional non-gravitational heating such as by Type la SNe. However, heating by Type la SNe is very likely to 
be important in this galaxy showing clear Fe enrichment and cool-hot gas interaction. Therefore, additional cooling 
processes must exist. 

The low temperature and high Fe abundance of the hot gas in the inner halo of NGG 5866 strongly suggest that 
the presence of cool gas in gas-poor galaxy bulges could trigger cool-hot gas interaction, resulting in a lower hot 
gas temperature and a larger fraction of SN ejecta detected in soft X-ray around the galaxy. However, most of the 
metal-enriched SNe ejecta is still missing from a clear detection in X-ray. 


Li Jiang-Tao would like to acknowledge Roberts Shawn R., Wang Q. Daniel, and the Suzaku helpdesk for discussions 
on the Suzaku background and also the anonymous referee for very helpful comments. This work is supported by NASA 
through the grant NNX13AE87G and NNX15AM93G. 
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